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Abstract. We report on the first three-dimensional hy- 
drodynamic simulations of secular and dynamical non- 
axisymmetric instabilities in collapsing, rapidly rotating 
stellar cores which extend well beyond core bounce. The 
resulting gravitational radiation has been calculated using 
the quadrupole approximation. 

We find that secular instabilities do not occur dur- 
ing the simulated time interval of several 10 ms. Models 
which become dynamically unstable during core collapse 
show a strong nonlinear growth of non-axisymmetric in- 
stabilities. Both random and coherent large scale initial 
perturbations eventually give rise to a dominant bar-like 
deformation (exp(±im0) with m = 2). In spite of the pro- 
nounced tri-axial deformation of certain parts of the core 
no considerable enhancement of the gravitational radia- 
tion is found. This is due to the fact that rapidly rotating 
cores re-expand after core bounce on a dynamical time 
scale before non-axisymmetric instabilities enter the non- 
linear regime. Hence, when the core becomes tri-axial, it 
is no longer very compact. 
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1. Introduction 

The rotation rate of a self-gravitating fluid body is usu- 
ally quantified by the parameter (3 :— E rot /\E pot \, where 
-E ro t is the rotational energy and E pot the potential en- 
ergy of the body (e.g., Tassoul 1978). It is well known from 
linear stability analysis that global non-axisymmetric in- 
stabilities can develop in a rotating body when j3 is suf- 
ficiently large (e.g., Chandrasekhar 1969; Tassoul 1978). 
Rotational instability arises from non-radial azimuthal (or 
toroidal) modes exp(±im0), where <j> is the azimuthal 
coordinate and where m is characterizing the mode of 
perturbation. The mode with m = 2 is known as the 
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"bar" mode. Driven by hydrodynamics and gravity, a ro- 
tating fluid body becomes dynamically unstable to non- 
axisymmetric perturbations on a time scale of approxi- 
mately one rotation period when (3 > /3d yn - Secular insta- 
bilities develop due to dissipative processes such as grav- 
itational radiation reaction (GRR; Chandrasekhar 1970), 
viscosity (Roberts & Stewartson 1963), or a combination 
of both (Lindblom & Detweiler 1977), and grow on a 
time scale of many rotation periods (Schutz 1989) when 

Pscc < (3 < 0dya- 

The values of the critical rotation parameters (3 sec and 
/3dyn depend on the density stratification and the angular 
momentum distribution of the rotating body. For Mac- 
Laurin spheroids, i.e., for incompressible, rigidly rotating 
equilibrium configurations, secular and dynamical bar in- 
stabilities set in for (3 SCC = 0.1375 and /3d yn = 0.2738, 
respectively (e. g., Chandrasekhar 1969; Tassoul 1978). In 
differentially rotating polytropes (p oc g 1+1 / n ) both the 
angular momentum distribution and to a lesser extent 
the polytropic index n affect the value of (3 SQC at which 
the secular bar mode sets in (Imamura et al. 1995). For 
rt = 3/2 polytropes and m — 2 perturbations Pickett et 
al. (1996) find that the classical dynamical stability limit 
/Sdyn ~ 0.27 holds, if the angular momentum distribution 
is similar to those of MacLaurin spheroids. However, con- 
siderably lower values of /3dyn are found when the latter 
gives rise to extended Keplerian-disk-like equilibria. 

It is unclear to what extent these results obtained an- 
alytically for MacLaurin spheroids and numerically for 
compressible configurations are applicable to collapsing 
cores, since they are derived for stationary equilibrium 
models. It has been shown, however, that during axisym- 
metric rotational core collapse conservation of angular mo- 
mentum can lead to very rapidly rotating configurations, 
whose rotation rates exceed f3 scc and even /3dyn (Tohlinc 
1984; Eriguchi & Miiller 1985; Monchmeyer et al. 1991; 
Zwerger & Miiller 1997). But whether such super-critical 
rotation in collapsing cores indeed leads to the growth of 
non-axisymmetric instabilities has not yet been demon- 
strated in hydrodynamic simulations. 

Besides of possibly being relevant for the collapse dy- 
namics of rapidly rotating stellar cores, tri-axial instabil- 
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ities have also been envisaged to boost the gravitational 
wave signal from rotational core collapse (for a review, see 
e.g., Piran 1990; Thornc 1995). This would be of great 
importance for the four long-baseline laser interferomet- 
ric gravitational wave detectors (GEO600, LIGO, TAMA, 
VIRGO) which will become operational within the next 
few years. They are designed to achieve sensitivities (for 
single bursts) down to gravitational wave amplitudes of 
\h\ <~ 10~ 22 (see e.g., Abramovici et al. 1992). In order to 
obtain a few events per year, one has to detect all core col- 
lapse supernovae out to the Virgo cluster of galaxies (dis- 
tance ~ 10 Mpc). However, simulations of axisymmetric 
rotational core collapse predict gravitational wave ampli- 
tudes of at most \h\ ~ 10~ 23 at that distance (Miiller 
1982; Finn & Evans 1990: Monchmeyer et al. 1991; Ya- 
mada & Sato 1995; Zwcrger & Miiller 1997). Significantly 
larger gravitational wave signals could be produced if non- 
axisymmetric instabilities due to rotation act in the col- 
lapsing core. According to this idea, a rapidly spinning 
core will experience a centrifugal hang-up and will be 
transformed into a bar-like configuration that spins end- 
over-end like an American football, if its rotation rate ex- 
ceeds the critical rotation rate(s). One has further specu- 
lated, that the core might even break up into two or more 
massive pieces, if (3 > /3dyn (e.g., Bonnell & Pringle 1995). 
It has been suggested that the resulting gravitational ra- 
diation could be almost as strong as that from coalescing 
neutron star binaries (Thorne 1995). The actual strength 
of the gravitational wave signal will sensitively depend 
on (i) the radius at which the centrifugal hang-up occurs 
and (ii) what fraction of the angular momentum of the 
non-axisymmetric core goes into gravitational waves, and 
what fraction into hydrodynamic waves. These sound and 
shock waves are produced as the bar or lumps, acting like a 
twirling-stick, plow through the surrounding mass layers. 

Recently, both semi-analytic and numerical methods 
have been used to compute the gravitational radiation 
produced by non-axisymmetric instabilities in rapidly ro- 
tating stars. Lai & Shapiro (1995) have studied secular 
instabilities in rapidly rotating neutron stars and have 
computed the resulting gravitational radiation using lin- 
earized dynamical equations and a compressible ellipsoid 
model. Houser et al. (1994), Houser & Centrella (1996) 
and Smith et al. (1996) have used three-dimensional hy- 
drodynamic codes to simulate the nonlinear growth of the 
dynamical tri-axial instability in rapidly (/3 = 0.3) rotat- 
ing polytropes (n = 3/2, n = l,n = 1/2). Scaling their 
results to neutron star dimensions (i.e., a polytrope with 
mass M <~ 1.4 M Q and radius R <~ 10 km) they have 
also calculated the gravitational radiation from the non- 
axisymmetric instability. They obtained a maximum di- 
mcnsionless gravitational wave amplitude \h\ <~ 2 10 -22 
for a source at a distance of 10 Mpc, the energy lost to 
gravitational radiation being AE <~ 10~ 3 Mc 2 . 

Concerning investigations, like those just discussed, it 
is important to note that they can be used to predict the 



gravitational radiation from rapidly rotating, stationary 
neutron stars, which might (or might not!) form as a con- 
sequence of rotational core collapse. They are not appro- 
priate, however, for predicting the gravitational wave sig- 
nature of the collapse of rapidly rotating stellar cores. We 
stress this difference here, because it is often overlooked. 

Up to now three-dimensional hydrodynamic simula- 
tions of rotational core collapse, which can follow the non- 
linear growth of non-axisymmetric instabilities during col- 
lapse, have only been performed by Bonazzola & Marck 
(1993, 1994), and by Marck & Bonazzola (1992). For their 
simulations they used a pseudo-spectral hydrodynamic - 
code and a polytropic equation of state. They computed 
the evolution of several initial models and found that the 
gravitational wave amplitude is within a factor of two of 
that of 2D simulations for the same initial deformation of 
the core (Bonazzola & Marck 1994). However, their sim- 
ulations were restricted to the pre-bounce phase of the 
collapse, and thus are less relevant. 

In the following we report on the first three-dimensio- 
nal hydrodynamic simulations of non-axisymmetric insta- 
bilities in collapsing, rapidly rotating stellar cores which 
extend well beyond core bounce. In addition to the dy- 
namics, we have also computed the gravitational radiation 
emitted during the evolution. Our study is a continuation 
and an extension of the recent work of Zwerger (1995) and 
Zwerger & Miiller (1997), who have performed a compre- 
hensive parameter study of axisymmetric rotational core 
collapse. The initial models, the equation of state and the 
hydrodynamic method are adopted from their study. From 
their set of 78 models we have only considered those which 
exceed the critical rotation rate(s) during core collapse. 

The paper is organized as follows: In section 2 we pre- 
sent the computational procedure used to follow the hy- 
drodynamic evolution of the core and to compute the grav- 
itational wave signal. Subsequently, in section 3, we first 
discuss the results of a two-dimensional simulation, which 
serves as a reference point for the three-dimensional runs. 
We then present various aspects of the evolution of the 
non-axisymmetric models and discuss the gravitational 
wave signature of the models. Finally, in section 4, we sum- 
marize our results, and conclude with a discussion of the 
shortcomings of our approach. 

2. Computational procedure 

2.1. Hydrodynamic method and equation of state 

The simulations were performed with the Newtonian mul- 
ti-dimensional finite- volume hydrodynamic code PROME- 
THEUS developed by Bruce Fryxell and Ewald Miiller 
(Fryxell et al. 1989). PROMETHEUS is a direct Eulerian 
implementation of the Piecewise Parabolic Method (PPM) 
of Colella & Woodward (1984). For the three-dimensional 
simulations we used a variant of PROMETHEUS due to 
Ruffert (1992), which utilizes multiple- nested refined equi- 
distant Cartesian grids to enhance the spatial resolution. 
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We do not include any general relativistic effects for the 
fluid and have neglected effects due to neutrino transport. 

Matter in the core is approximated by a perfect fluid 
using a simplified analytic equation of state (Janka et 
al. 1993; see also Zwerger & Muller 1997 for details of the 
implementation) . The pressure is a function of the density 
g and energy density u of the core matter and consists of 
two parts: 



P(g,u) = P P + P 



th ■ 



The polytropic part 
P P = K- Q r », 



(1) 



(2) 



which depends only on density, describes the pressure 
contribution of degenerate relativistic electrons or that 
due to the repulsive nuclear forces at high densities. The 
thermal part P t h mimics the thermal pressure of shock- 
heated matter. It is computed from the corresponding in- 
ternal energy density Uth by 



Pth = (Ah - l)uth, 



with 



r th = i.5. 



(3) 



The thermal energy density u t h in turn is given by the 
total energy density u through the relation 



(4) 



where u p is the energy density of the degenerate elec- 
tron gas. To account for the repulsive part of the nuclear 
forces the "polytropic" adiabatic exponent T p (see Eq. (||)) 
is increased from a value Ti close to 4/3 (resembling de- 
generate, relativistic electrons) to a value T2 — 2.5, if the 
density of a fluid element exceeds nuclear matter density 
which is set to g nuc — 2.0 10 14 gcm~ 3 . 

2.2. Gravitational wave emission 

In the transverse-traceless (TT) gauge a metric pertur- 
bation h can be decomposed as (e.g., Misner et al. 1973, 
chap. 35) 



/i x e> 



(5) 



with the unit linear polarization tensors defined, in 
spherical coordinates (r,9,<f), as 



e + = eg <X> e g - e^ 
e x = e B ®e,,+ e. 



(6) 



where eg and are the unit vectors in the correspond- 
ing coordinate directions and ® is the tensor product. 
The two independent polarizations h + and h x are calcu- 
lated from the fluid variables in Post-Newtonian approx- 
imation, keeping only the mass-quadrupole in the multi- 
pole expansion of the field. Let us consider an observer 



located at coordinates (r, 9, <fi) in a spherical coordinate 
system whose origin coincides with the center of mass of 
the core. Then the two independent polarizations are given 
by (e.y.,Zhuge et al. 1994) 



h x = 



CI/.. 

— - y ( I xx cos 2 <f> + t yy sin 2 4> + f xy sin 2(f)) cos 2 9 
+i zz sin 2 9 — (l xz cos <f> + f yz sin <fi) sin 29 



f XX Sin (f) - fyy COS <f> + l X y SVH 2(f) ) (7) 



G2 /l 

c 4 r \ 2 

+(I XZ smcj)— f yz cos0) sin# 



( 2^w ~ ixx) sin 2<^> cos 6 + f xy cos2</>cos( 



(8) 



Here the quantities are the second time derivatives 
of the Cartesian components of the reduced quadrupole 
moment tensor 1 . They are defined according to 



— / g(x,t-r/c) 



d 3 x, (9) 



where we have used the Einstein sum convention (i.e., 
summation over repeated indices) and where Sij is the 
Kronecker symbol. Note that the origins of the spherical 
coordinate system (used to specify the observer's position) 
and the Cartesian coordinate system (used to describe the 
matter distribution of the source) coincide. The coordinate 
systems are oriented relative to each other such that 9 = 
corresponds to the positive z-direction. 

We assume that the source possesses equatorial sym- 
metry (i.e., it is symmetric with respect to the transfor- 
mation z— > — z), and hence f xz — i yz = 0. Inspection 
of the remaining terms in Eqns. (0) and (§) then yields 
9 e {0, it/2} and 9 = as necessary conditions for extrema 
of the waveforms h + and h x (considered as functions of 
the observer's position), respectively. Thus, we calculate 
h x along the positive z-axis and h + in the equatorial plane 
(at an azimuthal position cf> — 0) to yield the maximum 
stresses^} We have checked for all models that we do not 
miss possible additional (short term) maxima of h + that 
could be present for an observer at some azimuthal posi- 
tion cj) j^z in the equatorial plane but could be overlooked 
by an observer located at <fr = 0. 

The quadrupole part of the total gravitational wave 
energy is given by 



E, 



GW 



1G 

5^ 



t ij i ij . 



(10) 



Instead of numerically approximating Eq. (|9|) with 
standard finite differences, it is superior (e.g., Finn & 



1 For all models discussed below, the value of h-\ 
turned out to be negligibly small. 
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Evans 1990; Monchmeyer 1993) to use an equivalent ex- 
pression derived independently by Nakamura & Oohara 
(1989) and by Blanchet et al. (1990): 



STF <2 g- (uV - x'd^) d 3 x 



(11) 



where the symmetric and trace free part of a doubly 
indexed quantity A^ is defined by 



STF{A i n := -A ij + -A jl - -8 ij A u 
1 1 2 2 3 



(12) 



The spatial derivatives dj$> of the gravitational poten- 
tial appearing in Eq. ( pi] ) are approximated by centered 
finite differences. 

For two reasons no back-reaction of the gravitational 
wave emission on the fluid (GRR) has been implemented. 
Firstly, the peak luminosity of gravitational radiation that 
can be expected from rotational core collapse is at most 
of the order of 10 50 ergs" 1 (Muller 1982; Finn & Evans 
1990; Monchmeyer et al. 1991; Zwerger & Muller 1997). 
When comparing this luminosity with the typical ener- 
gies (kinetic, rotational, potential, internal) E ~ 10 52 erg 
of the core one finds a GRR time scale of the order of 
100s, which is at least by a factor of 10 4 larger than the 
time interval over which we can calculate the core's evo- 
lution. Secondly, total energy conservation is violated by 
an amount AE ~ 10 49 erg in the best resolved 2D calcu- 
lation and by an amount AE ~ 10 50 erg in the 3D runs. 
This violation means an acceptable 1%-error as far as the 
dynamics is concerned. However, it dominates the change 
of energy radiated due to GRR (Eqw ~ 10 47 erg) in the 
time interval of our simulations by about a factor of 10 3 . 

2.3. Two-dimensional simulations 

Closely following Zwerger & Muller (1997), the axisym- 
metric (reference) simulations were performed in spheri- 
cal coordinates assuming equatorial symmetry. Thus, the 
angular grid covered the range 9 G [0,7r/2]. The number 
of equidistant angular zones was n§ = 18, 45 or 90, which 
corresponds to an angular resolution of 5°, 2° or 1°, re- 
spectively. The moving radial grid consisted of 360 zones. 
The grid was moved in such a way that the inner core 
(i.e., the subsonic inner part of the core) was always re- 
solved by 180 radial zones. The radial zoning varied from 
about 0.5 km in the unshocked inner core to about 50 km 
in the outermost parts of the outer core (for more details 
see Zwerger & Muller 1997). 

As in Zwerger & Muller (1997), the gravitational wave 
amplitude was computed by expanding the gravitational 
field into "pure-spin tensor harmonics" (Thorne 1980). 
The only non-vanishing quadrupole contribution is then 
given by 



A E2 _ 
^20 — 



G 16tt 3 / 2 d 2 



^_l") ( . i(1:(1/ . mj 



which is related to the (dimensionless) metric pertur- 
bations of Eq. (I5J) by 



1 /15 



^20 



h x =0 



(14) 



Using the counterpart of Eq. (11) in spherical coordi- 
nates (e.g., Zwerger & Muller 1997, Eq. 20), the second 
time derivative as well as the r 2 -weighting of mass ele- 
ments, both numerically troublesome (e.g., Finn & Evans 
1990), can be eliminated. 

In the quadrupole approximation the total gravita- 
tional wave energy is given by 



E, 



GW 



C A 1 

G 327 



A AE2 
a/i 20 

dt 



dt. 



(15) 



Higher-order terms in the multipole expansion of the 
gravitational wave field are negligible both for the wave 
amplitudes and the amount of radiated energy in the ax- 
isymmetric case (e.g., Monchmeyer et al. 1991, Zwerger 
1995). 

2.4- Three-dimensional simulations 

The three-dimensional calculations have been performed 
on multiple-nested refined equidistant Cartesian grids (see 
Ruffert (1992) for details of the method). As in the 2D 
calculations we assume the equatorial plane z — to be 
a plane of symmetry. Six centered grids were used each of 
which was refined by a factor of two. In direction of the 
rotation (z-) axis each grid had 32 zones, while the perpen- 
dicular grid planes consisted of 64 x 64 zones. The coarsest 
grid had a spatial extent of 2560 km covering the entire 
core with cubes of edge size 40 km. The finest grid cov- 
ered the innermost 80 km with a linear resolution of about 
1 km. Thus, the resulting spatial resolution was compara- 
ble to that of the 2D simulations. 

For the three-dimensional simulations we considered 
axisymmetric models of Zwerger & Muller (1997), which 
are stable against non-axisymmetric perturbations in the 
early stages of collapse (i.e., j3(t) < 0.1), but later exceed 
the critical rotation rate(s). Test calculations showed that 
non-axisymmetric instabilities did not grow in perturbed 
axisymmetric models with [3 ^5 0.1 on time scales of a 
few 10 ms. Using trilinear interpolation an axisymmetric 
model was mapped onto the three-dimensional grid when 
its rotation parameter reached a value of ~ 0.1. Typically, 
this occured a few milliseconds before core bounce. 

We then imposed certain non-axisymmetric perturba- 
tions and followed the subsequent dynamical evolution of 
the core in three spatial dimensions. Due to the cubic cells 
of our computational grid, the mapping procedure itself 
introduces deviations from axisymmetry that are symmet- 
ric with respect to the coordinate transformations 



0-> + 71-/2 



and 



(16) 
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Fig. 1. Snapshots of the density distribution (in units of [gem -3 ]) and the velocity field in a meridional plane for 
the axisymmetric model A4B5G5. The z-axis coincides with the axis of rotation (and symmetry). The contours are 
logarithmically spaced with intervals of 0.25 dex, they are shaded with darker grey values for higher density regions 
and labeled with their respective values. The time of the snapshot and the velocity scale are given in the top right 
corner of each panel. 



This implies that only azimuthal modes exp(im<^) with 
even values of m are contained in the "grid-mapping 
noise" . The amplitudes arc of the order of a few percent 
measured as relative deviations Sip := (tp — (ip)^) / (ip)<p of 
a scalar field ip(r,0,<j)) from its azimuthal mean (tp)^ := 
(27r) _1 J Q 27r ip(r, 9, 4>) d(j). This defines a lower limit for the 
size of the perturbation amplitudes, which we can impose 
explicitly onto the axisymmetric models. We point out 
that density perturbations 5g ~ 1% are about two orders 
of magnitude larger than those typically used in numerical 



stability analysis (e.g., Houser et al. 1994; Aksenov 1996; 
Smith et al. 1996). 

We used explicit random perturbations of the density 
with an amplitude of 10%, i.e., the density distribution of 
the mapped model was modified according to 

g(r, 9, </>) -> g(r, 9, <j>) • [l + 0.1 f(r, 9, 4>)] , (17) 

where / denotes a random function distributed uni- 
formly in the interval [—1,1]. All azimuthal modes that 
can be resolved on a given grid are therefore present in 
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the initial data, the m — 4 and m — 2 "grid-modes" 
mentioned above presumably being dominant. The cor- 
responding model is refcred to in the following as model 
MD1. 

In a second model MD2, we imposed in addition a large 
scale m = 3 azimuthal perturbation with an amplitude of 
5%, i.e., 

g{r, 6, 4>) -> g(r, 9, <j>) ■ [l + 0.05 sin(0.6 + 30)] . (18) 

Note that for convenience of notation we have given 
the form of the perturbation in spherical coordinates, al- 
though it is introduced on a Cartesian grid. 

3. Results 

Zwerger & Miiller (1997) have studied the overall dynam- 
ics of axisymmetric rotational core collapse and calcu- 
lated the resulting gravitational wave signal. In partic- 
ular, they have investigated how the dynamics and the 
gravitational radiation depend on the (unknown) initial 
amount and distribution of the core's angular momentum, 
and on the equation of state. Among the 78 models com- 
puted by them, Zwerger & Miiller (1997) found only one 
model (A4B5G5), whose rotation rate parameter consid- 
erably exceeds /?dyn near core bounce. However, /? remains 
larger than /?d yn for only roughly one millisecond, because 
the core rapidly re-expands after bounce and hence slows 
down. In addition, Zwerger & Miiller (1997) found three 
more models (A3B5G4, A3B5G5 and A4B5G4) that ful- 
filled (3 > p scc for several 10 ms (see Fig. 4 of Zwerger & 
Miiller 1997). 

The two "classical" dissipation mechanisms that can 
in principle drive secular non-axisymmetric instabilities, 
namely viscosity of core matter (Roberts & Stewartson 
1963) and gravitational radiation reaction (GRR; Chan- 
drasekhar 1970), are negligible during the time scales we 
have simulated evolution of the core (a few 10 ms). Nev- 
ertheless, some inner fraction of the collapsing core could 
still become unstable, because of acoustic coupling with 
or due to advection of matter into the outer core (see 
e.<7. ; Schutz 1983 and references therein). Therefore, we 
have experimented with two of the three models exceed- 
ing f3 sec (but not Pdyn) for several 10 ms by imposing non- 
axisymmetric perturbations on different hydrodynamic 
quantities (g, P, v) of different spatial character (m = 2, 
random) and amplitude (5%, 10%) at different epochs of 
the evolution. On a time scale of several 10 ms we have 
found neither indications for secular instabilities nor a 
significant enhancement of the gravitational wave signal 
(compared with the corresponding axisymmetric model). 

Thus, the only model remaining of the set of models of 
Zwerger & Miiller (1997), which is a promising candidate 
for the growth of tri-axial instabilities during its gravita- 
tional collapse, is model A4B5G5 (see above). Its prop- 
erties are described in the next subsection. Some global 
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t [10~ 3 s] 

Fig. 2. Radii of selected mass shells (for their defini- 
tion see Eq. ( po|) in the text) as a function of time for 
two three dimensional models MD1 (dashed lines) and 
MD2 (dashed-dotted lines), and the axisymmetric model 
A4B5G5 (solid lines). The mass inside the shown mass 
shells ranges from 0.1 M© to 0.9 M© in intervals of 0.2 M©. 
Lines corresponding to 0.1 M© and 0.7 M© are drawn in 
bold. 

quantities of the model have already been published by 
Zwerger & Miiller (1997). 

3.1. The axisymmetric model A4-B5G5 

All pre-collapse models of Zwerger & Miiller (1997) are 
axisymmetric n — 3 polytropes (p oc g 1+1 / n ) in rotational 
equilibrium whose angular velocity Q only depends on the 
distance w from the axis of rotation according to the ro- 
tation law 

QH K (l + (J 2 )" 1 . (19) 

Collapse is initiated by reducing the adiabatic index T 
from the value 4/3, which is consistent with the structure 
of a n = 3 polytrope, to a somewhat smaller value T\ (see 
Eq- §))■ 

The initial model A4B5G5gis the most extreme one in 
the large parameter set considered by Zwerger & Miiller 
(1997). It is the fastest (B\ = 0.04) and most differentially 
(A = 10 7 cm; see Eq. (pj|)) rotating model being evolved 
with the softest equation of state (T\ — 1.28). According to 
the criterion derived by Ledoux (1945; Eq. 77) for rotating 

2 If not stated otherwise, quoted numbers always refer to 
the 2D simulation with the best angular resolution, i.e., the 
one performed with ne = 90 (see section 2.3). 
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Fig. 3. Rotation rate parameter (3 as a function of time 
for the three-dimensional models MD1 (dashed) and MD2 
(dashed-dotted) . For comparison, the evolution of /3 of the 
two-dimensional model is shown, too (solid). The solid 
horizontal lines mark the critical values (3 SCC (lower line) 
and /?dyn (upper line), where MacLaurin spheroids become 
secularly and dynamically unstable to non-axisymmetric 
bar-like perturbations. 
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(Fig.[l]d). Like the outer shock, the inner shock eventually 
becomes prolate, too. 
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Fig. 4. Maximum density on the grid as a function of 
time for the three-dimensional models MD1 (dashed) and 
MD2 (dashed-dotted). For comparison we also show the 
maximum density of the two-dimensional model A4B5G5 
(solid). The solid horizontal line marks the assumed nu- 
clear matter density g n uc = 2 10 gem" 3 . 



stars, a rotation rate of (3 = 0.04 requires T\ > 1.3 for the 
star to be stable against the fundamental radial mode, 
i.e., collapse. The resulting model has a mass of 1.66M , 
a total angular momentum of L z = 3.85 1 49 ergs and an 
initial equatorial radius of 1280 km. 

Fig. 0a shows that model A4B5G5 has a torus- like den- 
sity stratification with an off-center density maximum, 
which is located in the equatorial plane 200 km away from 
the rotation axis. During collapse a kind of "infall chan- 
nel" forms along the rotation axis, a structure which has 
already been found in some models of Monchmeyer et 
al. (1991). There matter falls almost at the speed of free- 
fall towards the center. With decreasing latitude (i.e., in- 
creasing polar angle 0) matter is increasingly decelerated 
by centrifugal forces (Fig. 0b). At t m 28 ms, which is 
about 3 ms before core bounce, an oblate shock is visible 
(Fig. 0c). It begins already to form at the bottom of the 
infall channel 7 ms before core bounce (Fig. 0b). Pressure 
gradients decreasing with increasing polar angle lead to 
an asymmetric propagation of the shock its speed being 
larger in polar than in equatorial regions. The result is a 
"prolate shock in an oblate star" (Finn & Evans 1990). 

A second shock, which forms at the edge of the torus 
shortly before the density in the core reaches its maximum 
value inside the torus, propagates behind the outer shock 



The overall dynamics of the core can be conveniently 
analyzed in terms of the radial extent of different "mass 
shells" M(R) defined by 

M(R):=2ir / (g(r, 0)) r 2 sin0d0dr , (20) 
Jo Jo 

where ()^ denotes the azimuthal mean of the density 
to be calculated in the 3D simulations. M(R) is the mass 
contained in a sphere of radius R. Note that in more than 
one spatial dimension these mass shells are not identical 
with the well known Lagrangian mass shells. 

Inverting M(R) one obtains the radius Rm of the mass 
shell. This quantity is plotted as a function of time for dif- 
ferent masses M in Fig.0. One recognizes that maximum 
compression is reached in model A4B5G5 at t = 30.67 ms, 
which is followed by an expansion lasting for about 2 ms. 
Without any further significant oscillations the inner core 
evolves towards its new equilibrium: The radii of mass 
shells with M 0.7 Mq are nearly constant with time 
for t 33 ms (Fig. ||) . The same behavior can be deduced 
from the evolution of the rotation parameter (j3 tx for 
a homogeneous sphere with radius R) and of the maximum 
density on the grid (Figs.|| and 0). Figure0 also shows 
that the central core M(R) < 0.7 Mq oscillates with a 
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single dominant volume mode, which is in agreement with 
the results obtained by Monchmeyer et al. (1991). They 
found single volume modes in cores with coherent motion 
in equatorial and polar directions. This certainly occurs 
in model A4B5G5, since matter in the torus which con- 
tains most of the mass shows almost no motion in polar 
directions at all. 

Matter in the torus (about 0.7 Mq) stays in sonic 
contact, i.e., it contracts coherently. It also remains un- 
shocked well beyond core bounce. Hence, the torus resem- 
bles the qualitative features of a homologously contract- 
ing inner core that has been found analytically for the 
spherical collapse of polytropic cores (Goldreich & Weber 
1980, Yahil 1983) and numerically in slowly rotating cores 
with ellipsoidal density stratification (Finn & Evans 1990; 
Monchmeyer et aL1991). 

The main features of the quadrupole amplitude A^q 
calculated for model A4B5G5 are: A slow monotonic in- 
crease with time for the first 25 ms is followed by a pro- 
nounced negative spike at the time of bounce with no sub- 
sequent short-period oscillations (see Fig.^). We therefore 
consider the gravitational waveform to be closest resem- 
bled by signal type II, which is characterized by prominent 
spikes arising at the bounce(s) due to single dominant vol- 
ume mode(s) (Monchmeyer et al. 1991). 
The non-vanishing polarization, calculated for a source 
at lOMpc and for an observer located in the equa- 
torial plane of the core, has a peak value of h + = 
—3.5 10~ 23 . The energy of the quadrupole radiation is 
Egw = 7.8 10~ 8 Mq c 2 . Most of the spectral power is radi- 
ated at a frequency of v k, 200 Hz. This is the largest signal 
strength obtained for any axisymmetric collapse models in 
the set of Zwerger & Miiller (1997). However, the ampli- 
tude is still by at least one order of magnitude too small 
to be detectable even with the "advanced" LIGO inter- 
ferometer (see Abramovici et al. 1992 for the sensitivity 
limits), if the source is located in the Virgo cluster. 

3.2. Evolution of non-axisymmetric models 

The initial non-axisymmetric perturbations were intro- 
duced at t a = 28.13 ms (Figs.|]a and fjja). This is about 
2.5 ms before bounce (when the density reaches a maxi- 
mum inside the torus), and about 2ms before the core's 
rotation rates /3 *S 0.3, i.e., when it should become dynam- 
ically unstable. The additional perturbation of 5% ampli- 
tude and m = 3 azimuthal dependence imposed on model 
MD2 can hardly be distinguished from the 10% random 
noise perturbations of both models. The inner torus has 
contracted from an initial radius of 200 km to one of 60 km 
at t = to- Note that we refer to its density maximum, when 
we give its radial position. 

The subsequent rapid contraction of the rotating core 
is reflected by a steep rise of the maximum density, which 
peaks at t = 30.68 ms in both 3D models and in the ax- 
isymmetric one (Fig. |]) . At approximately the same time 



the contraction (Fig. the rotation parameter (3 (Fig. ||), 
the dissipation-rate of kinetic-infall energy into thermal 
energy and the gravitational potential energy reach their 
peak values. Thus, we consider t = 30.68 ms to be the time 
of bounce, although the central density does not reach its 
first peak then, which is usually considered as the bounce 
criterion. However, because of the torus-like density strat- 
ification of the models, the central density has less mean- 
ing. During the simulations it remained always much lower 
than the maximum density inside the torus (by 1-3 orders 
of magnitude). One also has to be careful, in particular in 
3D simulations, when deriving implications from "local" 
quantities. For example, one might be tempted to conclude 
from the fact that the maximum density exceeds nuclear 
matter density in model MD2, i.e.,g max > 210 14 gcm~ 3 
(Fig.|J), that this model contrary to model MD1 suffered 
a bounce due to the stiffening of the equation of state. 
However, nuclear matter density is exceeded only in a few 
zones inside the torus (see Fig. 0c), while the correspond- 
ing azimuthal average, which is determining the overall 
dynamics, is practically identical to that in model MD1. 
Therefore, we always compare the value of the maximum 
density with the corresponding azimuthal mean, before 
drawing any conclusions. 

Up to t w 33ms, i.e., up to about 2ms after 
bounce, the non-axisymmetric initial perturbations have 
not grown far enough to produce any significant changes 
in the overall dynamics of the collapsing core. This can be 
seen by comparing the time evolution of the radii of mass 
shells defined by Eq. ( p0| ) (see Fig.||). Also locally, non- 
axisymmetries are still small perturbations on a rapidly 
changing axisymmetric background. The density contrast 
in the torus is 5g 0.15 for model MD2 at t < 33 ms 
(Figs.^a-c). Hence, the distribution of hydrodynamic vari- 
ables in planes perpendicular to the rotation axis is still 
given by Fig.[j]. 

The subsequent evolution differs from that of the ax- 
isymmetric model and is also different between the two 
3D models. 

Model MD1: Due to the Cartesian geometry of the 
computational grid, the toroidal mode with m = 4 sticks 
out of the noise first. The inner torus assumes a transient 
square-like structure, which exists for « 7 ms. This does 
not influence the overall dynamics of the model signifi- 
cantly compared to the axisymmetric one, although we 
find a somewhat less compact density stratification when 
comparing the radii of the innermost mass shells with 
those of model A4B5G5 for t ^ 33 ms (Fig.|). When we 
stopped the simulation at t = 45.04 ms, the final configu- 
ration shows a prominent off-center density maximum in 
the equatorial plane rotating with a period of T « 5 ms at 
a distance of vj w 30 km. 

Model MD2: The development of the instability is 
shown in Figs.|5| and || By the time of bounce three dis- 
tinct density maxima are visible the density contrast being 
Sg ^5 0.15 inside the torus. These density maxima eventu- 
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Fig. 5. Snapshots of the density distribution (in units of [gem -3 ]) in the equatorial plane of model MD1. The contours 
are logarithmically spaced with intervals of 0.25 dex, they are shaded with darker grey values for higher density regions 
and labeled with their respective values. The time of the snapshot is given in the upper right corner of each panel. 



ally grow into three distinct clumps (Fig. 0c). The further 
evolution is characterized by intensive hydrodynamic ac- 
tivity produced by the twirling-stick action of the three 
clumps. Most notably are trailing "spiral arms" causing 
mass and angular momentum to be transported away 
from the non-axisymmetric, high density regions (Figs. 
|d, e). This is mirrored in the evolution of the innermost 
(M < 0.7 M Q ) mass shells, whose radii decrease during 
the time when the spiral arms are present. The largest de- 
viations from the monotonic expansion of the mass shells 
of model MD1 occur between t 33 ms and t ss 38 ms, 
when the spiral arms are most prominent in model MD2. 
Eventually, the three spiral arms merge (Fig. ^second last 



snapshot and Fig.fjje) and form a bar- like object inside the 
M = 0.7 M mass shell (Fig.| last snapshot and Fig.^f). 

3.3. Growth region of the instability 

For both 3D models we find that although perturbations 
have been imposed on the whole computational grid with 
the same relative amplitude initially, they grow only in 
the immediate vicinity of the inner torus, i.e., inside a 
cube with an edge length of w 100 km (Figs.|| and ^). 
This can be seen more clearly in Fig.pl where we have 
plotted the density (upper "curves" ) of each grid cell in 
the equatorial plane versus its (cylindrical) coordinate dis- 
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Fig. 7. Snapshots of the density distribution (in units of [g cm 3 ]) in the equatorial plane of model MD2. The contours 
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Fig. 6. Surfaces of constant density g — 10 12 gcm~ 3 for 
model MD2. The edges of the cubic box have a length 
of 160 km. The snapshots are taken at the same times as 
those in Fig. ^[ 

tance vj := y/ x 2 + y 2 from the origin. One notices consid- 
erable deviations from axisymmetry (which is mirrored in 
Fig. U by a large spread in density at a given distance m) 
only for w 40 km. That this is not an artefact due to 
the usage of a multiple-nested refined grid can be seen as 
follows: The cube-shaped boundaries of the single grids 
are located at a Cartesian coordinate distance Abound G 
{40 km, 80 km, 160 km, . . .} from the origin (£ £ {x, y, z}). 



Fig. 8. Scatter plots of the density (upper "curves", loga- 
rithmic scale, in units of [gem -3 ]) and the Mach number 
Ma* (lower "curves", calculated in a coordinate system 
rigidly rotating with the torus) in the equatorial plane 
z — versus distance vj := \/ x 2 + y 2 from the rotation 
axis at different instants in time (see upper right corner 
of each panel). The bold horizontal line separates super- 
sonic (upper) from subsonic (lower) flow regions. The thin 
vertical lines indicate the boundaries Abound of the nested 
grids (see text). 



Depending on the azimuthal angle 4> this implies for each 
grid a doubling of the (linear) zone size at cylindrical co- 
ordinate distances Abound G [Abound, (1 + \/2) • Abound]- If 

the growth of the perturbations was damped notably by 
the coarser numerical resolution at larger radii one would 
expect the spread in density as function of vj to change 
significantly at the limiting values of tUbound) which is ob- 
viously not the case in Fig.||. 

In contrast, we observe a striking correlation between 
the Mach-number Ma* calculated in a rotating coordinate 
system whose angular velocity fi w 1800s -1 is that of the 
torus at t = 30 ms (when (3 = /3dyn) and the domain where 
the density varies considerably with azimuthal angle: Non- 
axisymmetric perturbations grow only in regions where 
the flow is subsonic, i.e., where Ma* < 1 (Fig.||). 

Our qualitative interpretation of this result is as fol- 
lows: Non-axisymmetric instabilities have been shown to 
occur in MacLaurin spheroids as well as in a large va- 
riety of (differentially) rotating compressible fluid bodies 
in equilibrium, if the rate of rotation is sufficiently large. 
However, being equilibrium configurations, the instabil- 
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ities in these models grow in the absence of any radial 
motion. When viewed from a suitably chosen rotating co- 
ordinate system sonic contact along the azimuthal coordi- 
nate is established in these models and allows for a coher- 
ent growth of the instability. According to our results, the 
latter condition seems to be sufficient for the growth of 
non-axisymmetric instabilities also in collapsing rotators. 
If on the other hand sonic contact is not given along the 
azimuthal coordinate due to large radial velocities — like 
in the outer core of our models — it is difficult to imagine 
how global tri-axial deformations can develop coherently. 
Thus, we suppose that sonic communication in azimuthal 
direction is also a necessary condition for the growth of 
non-axisymmetric instabilities. 

One might argue that these results depend on our par- 
ticular choice of the rotating coordinate system. However, 
establishing co- rotation at any distance 137 (i.e., chosing a 
coordinate system rotating with an angular velocity equal 
to that of the fluid at that distance) outside the domain 
of the inner torus does not allow to "transform away" su- 
personic velocities in the collapsing (outer) core, since the 
radial component dominates the angular component of the 
velocity field by a large amount. 

3.4- Gravitational wave signal 

Although models MD1 and MD2 show prominent devia- 
tions from axisymmetry, the gravitational wave signal is 
changed only marginally compared with the axisymmet- 
ric calculation. In particular, the maximum wave ampli- 
tudes are equal for the 3D and the 2D simulations within 
an accuracy of 2%. The waveforms and the energy radi- 
ated in form of gravitational waves are displayed in Fig.^J. 
The peak amplitude of h + is reached at about the time of 
bounce in both the 2D and the 3D models. Up to this point 
in the evolution the waveforms are identical, too. This is 
not surprising, since the initial perturbations of models 
MD1 and MD2 have not grown significantly until the time 
of bounce, which occurs only 0.5 ms after /3dyn is reached. 
As discussed above, the density contrast is 5g £5 0.15 in- 
side the inner torus before bounce. Accordingly, h x , which 
is of genuine non-axisymmetric origin nearly vanishes until 
bounce (Fig.||). 

The maximum amplitudes of h x are reached during 
the further evolution, when the spiral arms merge to form 
the final bar (model MD2) or when the transient m = 4 
structure evolves to some "bar-like" configuration (model 
MD1). These maxima are only of the order of 10% of the 
maximum value of \h + \. The amplitudes of the cross- and 
plus-polarizations finally become comparable, because the 
inner core (torus) approaches its new rotational equilib- 
rium and thus the time derivatives of the quadrupole mo- 
ments due to radial motion become steadily smaller. The 
periodicity seen in h x is due to a nearly solid-body type 
rotation of the high density regions of the core. 
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Fig. 9. Gravitational wave amplitudes rh + (x = r,y = 
0, z = 0), (thick lines), rh x (x = 0,y = 0,z = r), (thin 
lines) and radiated energy -Egw as functions of time for 
the three-dimensional models MD1 (dashed) and MD2 
(dashed-dotted) in comparison with the axisymmetric cal- 
culation (solid). The peak values \rh\ ~ 1000cm corre- 
spond to dimensionless amplitudes \h\ w 3 10~ 23 for a 
source at r = 10 Mpc. 



Comparing the 3D results with the 2D signal, one no- 
tices that the 3D waveforms do not reproduce the small 
additional local minimum in h + visible at t k 31ms for 
the axisymmetric model. Also the (absolute values of the) 
time derivatives of h+ are smaller for the following 2 ms. 
This leads to a smaller value of the amount of energy ra- 
diated in form of gravitational waves (Eq.[l0|; for model 
MD1 h + is roughly proportional to i xx and f yy during 
this time interval because of the near axisymmetry of the 
core at this time). We have analyzed carefully the unex- 
pected fact that compared to the axisymmetric models 
we obtain only 65% of the gravitational wave energy for 
the two models MD1 and MD2, which during this epoch 
are approximately, but not exactly axisymmetric. Several 
tests and comparisons (e.g., with poorer resolved 2D cal- 
culations) have shown that this is not a possibly unno- 
ticed 3D effect, but due to a somewhat lower "angular" 
resolution of the 3D simulation compared with the best re- 
solved 2D run (see Rampp (1997) for details). The lower 
resolution gives rise to a larger violation of total energy 
conservation, most of which occurs during bounce. Since 
all dynamical quantities of model MD1 - local as well 
as global ones - agree with the 2D results within an ac- 
curacy of a few percent, we consider the observationally 
relevant waveforms calculated in the 3D simulations to be 
reliable within the underlying approximations (cf. sect. 2). 
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We also point out the fact that the squared time deriva- 
tives of the moments f jj enter the formula for the radiated 
energy. Small deviations in the waveforms therefore can 
account for quite large differences in the energies. It is not 
surprising, though not easy to explain in detail, that the 
3D waveforms differ from the axisymmetric waveforms at 
later times t <; 33 ms, when the dynamical evolution of 
these models is changed significantly due to their non- 
axisymmetric inner cores. 

Given the dynamical evolution of the core, the small 
magnitude of the amplitudes h x can easily be explained 
by utilizing the well known order-of-magnitude argument 
for the quadrupole waveforms 



rh y 



G_I_ 
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G Q.IMR 2 
~? T 2 



100 cm. 



(21) 



According to the results of our computations, we 
have inserted the quadrupole moment f of a homoge- 
neous bar with mass M ~ 0.5 M©, length R ~ 100 km 
and rotation period T ~ 5 ms. Approximating the non- 
axisymmetrically distributed mass of ~ 0.5 Mq as two 
point masses orbiting each other at a distance of 100 km 
with an angular velocity equal to the observed one, yields 
the same order-of-magnitude for h x . 

Finally, by using the quadrupole formulae (Eqns.^ and 
|8|) we may have underestimated the amount of gravita- 
tional radiation produced by model MD2, particularly 
during the phase when the m — 3 symmetry of the in- 
ner core is most prominent. We estimate the order-of- 
magnitude of the mass-octupole contribution to the signal 
by (e.g., Blanchet et al. 1990, Eq. 6.8) 
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(22) 



where v ~ 0.1 c and M ~ 0.5 M Q and a Schwarzschild 
radius R s ~ 1.510 5 cm was assumed. Hence, the mass- 
octupole radiation cannot account for significant enhance- 
ment of the peak amplitudes of our model MD2 compared 
with the axisymmetric case (Fig.^) or MDl, although it 
might change the details of the waveforms. 

3.5. Second collapse of a non-axisymmetric core 

After bounce, the proto-neutron star settles into its new 
rotational equilibrium, and cools and contracts on a sec- 
ular time scale, which is given by the neutrino-loss time 
scale (~ 10 s; see e.<?. ; Keil & Janka 1995 and references 
cited therein). According to linear stability analysis tri- 
axial perturbations grow on this time scale, provided the 
rotation parameter (3 0.14. For numerical reasons (be- 
cause our hydrodynamic code is explicit and thus the 
time step is limited by the CFL stability condition; see 
e.g.,LeVeque 1992) as well as for physical reasons (ne- 
glect of weak interactions and neutrino transport) our ap- 
proach is inadequate for simulating this secular evolution. 
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Fig. 10. Rotation parameter (3 as a function of time for 
the three-dimensional model MD3 (dotted), which suffers 
a second collapse. For comparison the evolution of (3 of 
model MDl is shown, too (dashed). 



However, we are able to consider a sudden reduction of 
the stabilizing pressure in a deformed, non-axisymmetric, 
rapidly rotating post-bounce core (model MD3). This ex- 
treme case cannot provide the answer to the question 
whether secular instabilities will indeed grow to nonlin- 
ear amplitudes and what their influence will be on the 
evolution of the core. Nevertheless, such a simulation can 
shed some light on the problem how much gravitational 
radiation can be expected, when a rapidly rotating non- 
axisymmetric neutron star forms. 

We choose model MDl at t = 43.12 ms as the starting 
point for the simulation (Fig.|i~l|a). At this time we (sud- 
denly) reduce Ti from its original value of 1.28 to a value 
of 1.2. This renders the rotating core unstable against a 
second dynamic collapse. In passing we note that the sta- 
bility criterion due to Ledoux (1945; Eq. 77) would require 
Ti 1.18 for a rigidly rotating core in equilibrium with 
(3 w 0.16 to collapse. 

The overall contraction, bounce and re-expansion of 
model MD3 is reflected in the time evolution of the rota- 
tion parameter (3. Note, that concerning the development 
of non-axisymmetric instabilities the absolute value of (3 is 
not very relevant here, because the initial "perturbations" 
are already in the non-linear regime. Compared to the 
(first) collapse of model MDl we find a somewhat longer 
time scale for the contraction and re-expansion in model 
MD3. The peak value of j3 m 0.23 is considerably smaller 
than that of model MDl (/3 0.34), because only mass 
shells with M ^ 0.5 M© contract to radii similar to those 
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Fig. 11. Snapshots of the density distribution (in units of [gem -3 ]) and the velocity field in the equatorial plane of 
model MD3. The contours are logarithmically spaced with intervals of 0.25 dex, they are shaded with darker grey 
values for higher density regions and labeled with their respective values. The time of the snapshot and the velocity 
scale are given in the upper right corner of each panel. 



reached during the collapse of model MD1. Mass shells of 
model MD3 with M ^ 0.8 Mq remain nearly unaffected 
by the sudden pressure reduction at t = 43.12 ms. 

During the collapse of model MD3 several pronounced 
density maxima are formed. They are all located inside 
the torus reaching peak values of g = 5.3 10 14 gcm -3 . The 
five initial clumps forming at the density maxima contain 
a mass of w O.1M each (Fig. |Tl|b) . During the further 
evolution their number decreases as they merge one after 
the other. Eventually, just two of them remain, which form 



a bar-like or binary-like deformed central region (Fig. |Tl|c 
and d). 

The maximum density on the grid (reached in the 
clumps) sharply rises from 2 10 13 gcm -3 to 5.3 10 14 gcm~ 3 
within a fraction of a millisecond after pressure reduction. 
The collapse is, however, rapidly decelerated by the action 
of centrifugal forces, which also eventually cause a bounce 
at t « 46 ms (Fig. [l0]). The pressure increase resulting from 
matter whose density exceeds nuclear matter density is 
dynamically unimportant, as only a negligible amount of 
mass inside some of the clumps is involved (Fig. [ll]) . 
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Fig. 12. Gravitational wave amplitudes rh + (x = r,y = 
0,z — 0), (thick lines) and rh y (x~Q,y = Q,z = r), (thin 
lines) as functions of time for the three-dimensional model 
MD3 (dotted), where r denotes the distance to the source. 
For comparison the amplitudes of model MD1 are shown, 
too (dashed). 



Model MD3 shows a highly non-axisymmetric, bar-like 
density stratification already before the (second) bounce 
occurs at t 46ms (Fig.[Tl"fc>). This is different from the 
behaviour of models MD1 and MD2, where large non- 
axisymmetries (quadrupole moments) only appear well af- 
ter bounce, i.e., after the most rapid phase of the evolu- 
tion. Accordingly, one expects larger time derivatives of 
larger quadrupole moments and therefore a stronger grav- 
itational wave signal from model MD3 than either from 
model MD1 or MD2. 

The actual waveforms for model MD3 are shown in 
Fig. |l2]. Compared to model MD1 one notices a consider- 
ably larger value for the genuinely non-axisymmetric po- 
larization amplitude h x . The maximum values \rh x \ 
800cm and \rh + \ « 300cm do, however, neither exceed 
the peak value of \rh+ \ obtained in the axisymmetric sim- 
ulation at bounce nor those of the three-dimensional mod- 
els MD1 and MD2 (Fig.[|). 

Additional (axisymmetric) simulations with Ti re- 
duced to 1.2 or 1.25 have shown that the peak value of 
/3, the time scale of the (second) collapse as well as the 
peak values of the gravitational waveforms do not depend 
strongly on the exact value of Ti . 



4. Summary and discussion 

We have presented three-dimensional hydrodynamic sim- 
ulations of the collapse of rapidly rotating stellar iron 
cores. The matter in the core has been described by a 
simple analytical equation of state. Effects due to neu- 
trino transport have been neglected. Hydrodynamics and 
self-gravity has been treated in Newtonian approximation. 
The initial models for our simulations have been taken 
from the comprehensive parameter study of axisymmetric 
rotational core collapse performed by Zwerger & Muller 
(1997). 

For our study we have selected those axisymmetric 
models which are secularly or dynamically unstable with 
respect to non-axisymmetric perturbations. Hence, in or- 
der to be selected the rotation rate parameter [3 of a 
model had to exceed the critical rotation rates /3 sec ~ 0.14 
or Pdyn ~ 0.27 during its evolution. These critical rota- 
tion rates have been derived for incompressible MacLaurin 
spheroids, but also hold approximately for a wide range 
of compressible equilibrium configurations (e.g., Tassoul 
1978). After mapping the rapidly rotating axisymmet- 
ric cores onto the three-dimensional computational grid, 
we imposed an initial non-axisymmetric perturbation and 
simulated the further evolution with a three-dimensional 
variant of the PROMETHEUS hydrodynamic code. The 
simulations cover a time interval from a few milliseconds 
before core bounce up to several tens of milliseconds after 
bounce. 

We have used the quadrupole formula to calculate the 
gravitational wave signal and have estimated the contribu- 
tion of the next-leading radiative-multipole order, which 
is negligibly small. Because of its smallness gravitational 
radiation reaction was not taken into account. 

Our 3D simulations show that in two models, which are 
secularly but not dynamically unstable (0.15 < [3 < 0.2) 
non-axisymmetric perturbations do not grow. Consistent 
with this result, we also observe no significant enhance- 
ment of the gravitational wave emission in these models 
compared to the axisymmetric case. 

Among the models investigated by Zwerger & Muller 
(1997) there is only one model where [3 > /3d yn during 
the evolution. This is their most rapidly and most differ- 
entially rotating model evolved with the softest equation 
of state. This axisymmetric model was perturbed about 
2 ms before bounce by imposing a random density dis- 
tribution of 10% amplitude. Despite the random initial 
perturbation, m = 2 and m = 4 toroidal modes are found 
to dominate the early evolution of the model, because of 
the set of cubic grids used in the simulation. Therefore, 
a second simulation was performed, where an additional 
tti — 3 toroidal perturbation of 5% amplitude was imposed 
initially. 

The gross features of the evolution are quite similar 
in the axisymmetric and non-axisymmetric simulations. 
In both cases the overall evolution is characterized by a 
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rapid contraction, bounce and re-expansion of the inner 
core. In the 3D simulations we, in addition, observe the 
growth of the initial perturbations. Until several millisec- 
onds (or dynamical time scales of the inner core) after 
bounce the density distribution resembles the symmetry 
of the initially dominating modes. Subsequently, the core 
evolves towards a bar-like shaped configuration. In none of 
the 3D models the gravitational wave amplitude exceeds 
the axisymmetric value of \h\ = 3.5 10 -23 (for a source 
at lOMpc), three dimensional effects on the waveforms 
being only of the order of 10%. Even when the collapsed, 
bar-like, rapidly rotating inner core is forced into a second 
collapse (by artificially reducing its adiabatic index), we 
do not observe gravitational wave amplitudes significantly 
larger than \h\ ~ 10~ 23 (for a source at lOMpc). 

In order to judge the implications of our results for ro- 
tational core collapse in general and for gravitational wave 
astronomy in particular, the following possible limitations 
of our approach should be kept in mind: 

(i) Since GM/Rc 2 ^ 0.2 for the axisymmetric models 
of Zwerger & Miiller (1997), general rclativistic ef- 
fects can be viewed as moderate corrections to New- 
tonian gravity as far as the collapse dynamics is con- 
cerned. However, since general relativity counteracts 
the stabilizing influence of rotation on radial modes 
(e.g., Tassoul 1978), its influence on the stability prop- 
erties of (even only moderately compact) rotating iron 
cores can be of considerable importance. When GR is 
taken into account, prc-collapsc models with a larger 
amount of angular momentum (and the same EOS) 
than in the Newtonian approximation can collapse. 
Furthermore, GR models will maintain higher densi- 
ties as well as larger rotation parameters for a longer 
time interval after bounce as compared to Newtonian 
models (see the discussion in Zwerger & Miiller 1997). 

(ii) The secular evolution (i.e., the evolution on time scales 
of hundreds of milliseconds) of the collapsed core is not 
known even in the axisymmetric case. Simulating this 
evolution requires models with detailed microphysics 
and thermodynamics, and an adequate treatment of 
the neutrino transport. However, even in 2D such sim- 
ulations are still prohibitively time consuming. Thus, 
there might exist additional initial models different 
from those presently available, which could become un- 
stable to tri-axial perturbations on a secular time scale 
(see however the next point). 

(iii) Concerning the limited set of axisymmetric initial mod- 
els, we adopt the arguments given in Zwerger & Miiller 
(1997). They claim that the large parameter space con- 
sidered in their study most probably comprises the 
whole domain where actual pre-collapse rotating iron 
cores are to be found. They further argue that using a 
realistic equation of state instead of a polytropic one 
will change the details of the evolution, but will not 
give rise to qualitatively different results. In particu- 



lar, it is unlikely that there exists a much larger set of 
pre-collapse initial models, which fulfill j3 > /3dyn and 
which are sufficiently compact for an extended time 
interval after core bounce. 

But even if there arc models which fulfill these con- 
ditions and become triaxial, we do not expect them 
to produce a considerably stronger gravitational wave 
signal than the "best" axisymmetric models. Provided 
the initial amount of angular momentum is not un- 
reasonably large, the critical rotation rates for non- 
axisymmetric instabilities can only be exceeded in col- 
lapsing cores with relatively soft equations of state 
(Eriguchi & Miiller 1985, Zwerger & Miiller 1997). 
But these soft equations of state drastically reduce 
both the mass and the radius of the inner core (Zwer- 
ger & Miiller 1997). However, according to our 3D 
calculations, it is only in the inner core, where non- 
axisymmetric perturbations can grow significantly. Ta- 
ken together these considerations suggest that bar-like 
inner cores have masses less than 0.5 M©, radii less 
than 100 km, and rotation periods greater than 1ms. 
This implies a maximum gravitational wave amplitude 
\h\ <, 10~ 22 (for a source at lOMpc). 

Given the most recent rates of core collapse supernova 
events (e.g., Cappellaro et al. 1997), which not all may in- 
volve a rapidly rotating core, we conclude that it is rather 
unlikely to expect gravitational wave signals from dynamic 
non-axisymmetric instabilities with sufficient strength and 
rate (i.e., several per year) to be observable with the large 
interferometric gravitational wave antennas presently un- 
der construction. 

Movies in MPEG format of the dynamical evolution of 
all models are available in the world- wide- web at http: / / 
www. mpa-garching. mpg. de/~ mjr / GRAV/grav3.h tml 
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